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ABSTRACT  '  V 

The  Optimum  Thermal  Interpolation  System  fOTIS  f  is  an  ocean  thermal  analysis  product  developed  for 
real-time  operational  use  at  the  ITS.  Navy  \s  f  leet  Numerical  Oceanography  Center!  It  functions  in  an  analysis- 
prediction-analysis  data  assimilation  cycle  with  an  ocean  mixed-layer  mode)  using  the  optimum  interpolation 
formulation  of  Alaka  and  Elvander  (  1972 1.  Thus,  climatology  serves  as  the  first-guess  field  for  thCanuKsis  of 
synoptic  ship,  bathythermograph,  buoy,  and  satellite  data,  with  the  prediction  from  the  mixed-layer  model 
treated  as  a  special  class  of  data. 

An  operational  test,  involving  comparisons  of  OTIS  against  an  existing  operational  ocean  thermal  structure 
model,  was  conducted  during  February.  March,  and  April  I98X.  Qualitative  comparison  of  the  two  products 
suggests  that  OTIS  gives  a  more  realistic  representation  of  subsurface  anomalies  and  horizontal  gradients.  Quan¬ 
titative  comparison  of  the  two  products  with  over  9900  independent  (i.e..  unassimilated)  bathythermograph 
observations  to  calculate  apparent  root-mean-squarc  <rms)  errors  indicate  that  OTIS  gives  a  more  accurate 
analysis  of  the  thermal  structure,  with  improvements  largest  below  the  base  of  the  mixed  layer.  Regional  differences 
in  the  relative  performance  of  the  models,  which  arc  probably  related  to  (he  validity  of  the  prescribed  statistical 
parameters  required  by  OT  IS,  are  also  noted  from  verification  against  the  bathythermograph  data.  OTIS  performs 
best  relative  to  the  existing  model  in  the  Indian  Ocean,  and  poorest  relative  to  this  mixiel  in  the  eastern  midlatitude 
Pacific,  in  the  latter  area,  however,  both  models  exhibit  their  lowest  error  levels  and  perform  best  relative  to 
climatology. 

'  Subgrid-scale  noise  contaminates  the  apparent  rms  error  statistics  and  obscures  the  relative  grid-scale  accuracy 
of  the  models.  A  general  technique  is  presented  to  estimate  grid-scale  rms  errors  which  arc  uncontaminaled  by 
subgrid-scale  and  instrumental  noise  in  the  observations.  Application  of  the  technique  indicates  that  the  grid- 
scale  errors  for  OTIS  are  typically  2(TV  less  than  those  of  the  other  model. 

OTIS  was  designated  as  the  Navy's  new  operational  global-scale  ocean  thermal  analysis  product  in  July  1988. 


I.  Introduction 

Real-time  monitoring  of  thermal  structure  in  the 
upper  ocean  is  becoming  increasingly  important  as  we 
seek  a  better  understanding  of  global  climate  change. 
In  addition,  an  accurate  representation  of  the  upper 
ocean  is  important  in  both  long-range  and  medium- 
range  weather  prediction,  fisheries  management,  and 
the  operation  of  underwater  acoustic  systems.  The  ll.S. 
Navy's  Fleet  Numerical  Oceanography  Center 
(FNOC).  Monterey.  California,  has  operated  com- 
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puterized  analyses  for  real-time  monitoring  of  ocean 
thermal  structure  on  global  and  regional  scales  for  over 
25  years.  In  general,  these  analyses  have  combined  real¬ 
time  ship,  bathythermograph  (bathy).  buoy,  and  sat¬ 
ellite  observations  with  climatology  and  predictions 
from  numerical  ocean  models  to  produce  fully  auto¬ 
mated  "nowcasts"  of  ocean  thermal  structure. 

The  Optimum  Thermal  Interpolation  System 
(OTIS)  is  the  latest  ocean  thermal  analysis  product 
developed  for  operational  use  at  FNOC.  OTIS  is  based 
on  the  Optimum  Interpolation  (OI)  data  assimilation 
technique  ofGandin  (  l%5).  which  is  used  widely  in 
meteorology.  The  OI  technique  has  been  applied  to 
oceanographic  analysis  for  research  applications  by 
Brcthcrton  ct  al.  (  1976  ).  Freeland  and  Gould  (  1976  I. 
McWilliams  (  1976).  White  (  1977).  White  and  Bern¬ 
stein  (1979).  Bernstein  and  White  (|9SI).  Clancy 
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(  1983),  Roemmich  (1983),  Robinson  and  Leslie 
(  1985),  McWilliams  et  al.  ( 1986),  Hua  et  al.  ( 1986), 
Robinson  et  al.  ( 1986),  Cummings  ( 1986),  Carter  and 
Robinson  ( 1987),  Robinson  et  al.  ( 1987).  and  others. 
Because  the  OI  technique  is  well  suited  for  handling 
data-sparse  areas  and  assimilating  data  with  differing 
error  characteristics,  it  is  an  attractive  approach  for 
operational  ocean  thermal  analysis. 

Basically  and  in  general,  the  Ol  technique  maps  ob¬ 
servations  distributed  nonuniformly  in  space  and  time 
to  a  uniformly  gridded  synoptic  representation,  or 
analysis,  of  the  target  field.  As  in  most  objective  analysis 
techniques,  the  concept  of  resolvable  versus  subgrid- 
scale  features  is  fundamental  in  OI.  Resolvable  features 
are  those  of  spatial  extent  greater  than  twice  the  mesh 
length  of  the  analysis  grid,  which  are  thus  capable  of 
being  at  least  marginally  represented  by  the  grid. 
Subgrid-scale  features  are  those  with  spatial  extent  less 
than  twice  the  mesh  length  which,  as  a  result,  cannot 
be  resolved  by  the  grid.  As  far  as  the  analysis  is  con¬ 
cerned,  subgrid-scale  features  contribute  to  the  error 
in  the  observations. 

An  OI  analysis  is  constructed  as  a  background  or 
first-guess  field  plus  anomalies  relative  to  this  field.  In 
general,  the  anomaly  at  a  particular  gridpoint  is  given 
by  a  weighted  combination  of  observed  anomalies,  with 
the  space-time  autocorrelation  function  for  the  re¬ 
solvable  anomalies  governing  w  hich  observations  con¬ 
tribute.  The  OI  technique  provides  the  optimum 
weights  applied  to  the  observed  anomalies  such  that 
the  resulting  analysis  error  will  be  minimized  in  a  least- 
squares  sense.  The  technique  also  provides  an  estimate 
of  this  error.  The  basic  inputs  to  this  process  are  the 
statistics  defining  both  the  resolvable  and  subgrid-scale 
variability  of  the  target  field  about  the  background  field 
and  the  instrumental  error  characteristics  of  the  mea¬ 
surement  system  providing  the  observations. 

In  operational  meteorological  applications  of  the  OI 
technique,  the  first-guess  field  for  the  analysis  is  given 
generally  by  a  model  prediction  from  the  previous 
analysis.  Climatology,  the  long-term  mean  state  of  the 
atmosphere,  does  not  enter  the  process.  For  operational 
oceanographic  applications,  however,  this  is  not  a  vi¬ 
able  approach:  the  present  generation  of  ocean  thermal 
prediction  models,  the  atmospheric  forcing  used  to 
drive  them,  and  the  observational  data  base  available 
to  update  them  are  simply  not  good  enough,  at  least 
for  global-scale  application,  to  keep  the  evolv  ing  ther¬ 
mal  fields  "on  track."  Climatological  ocean  thermal 
structure,  derived  from  historical  data,  must  be  used 
as  a  constraining  factor.  Consequently.  OHS  follows 
the  approach  of  Alaka  and  Flvandcr  (  1972  )  by  using 
climatology  as  the  first-guess  field  for  the  OI  analysis 
and  treating  input  from  an  ocean  prediction  model  as 
a  special  class  of  data.  The  prediction  model  used  by 
OTIS  is  the  Therrmxlynamic  Ocean  Prediction  System 
(TOPS)  mixcd-laver  model  described  bv  Olancv  and 
Poliak  (  1983). 


The  Expanded  Ocean  Thermal  Structure  (EOTS) 
analysis  has  been  operational  at  FNOC  for  almost  ten 
years  (Clancy  1987).  EOTS  combines  real-time  ocean 
thermal  observations  with  climatology  via  the  Fields- 
by-lnformation-Blending  (FIB)  analysis  technique 
( Holl  et  al.  1 979 )  to  produce  a  three-dimensional  anal¬ 
ysis  or  “nowcast”  of  ocean  thermal  structure.  A  number 
of  operational  EOTS  runs  are  made  daily  at  FNOC 
corresponding  to  various  regions  of  high  Navy  interest. 
When  coupled  additionally  in  certain  regions  to  the 
TOPS  mixed-layer  model  in  an  analysis-prediction- 
analysis  data  assimilation  cycle,  the  system  is  referred 
to  as  the  TOPS-Coupled  EOTS  (TEOTS)  analysis 
(Clancy  and  Poliak  1983:  Martin  et  al.  1985).  Thermal 
fields  from  EOTS  and  TEOTS  are  used  to  support 
FNOC  acoustic  predictions  and  transmitted  to  nu¬ 
merous  users  through  a  variety  of  communications 
links. 

An  operational  test  was  conducted  between  1  Feb¬ 
ruary  and  30  April  1988  to  establish  OTIS  as  a  valid 
replacement  for  the  global-scale  implementation  of 
TEOTS.  This  involved  qualitative  comparison  of  ther¬ 
mal  fields  produced  by  OTIS  with  those  produced  by 
TEOTS.  and  quantitative  verification  of  both  models 
against  independent  bathy  data.  During  the  test.  OTIS 
and  TEOTS  functioned  on  the  same  horizontal  grids 
(the  standard  FNOC  Northern  Hemisphere  and 
Southern  Hemisphere  63  X  63  polar  stereographic 
grids:  see  Clancy  and  Poliak  1983).  used  the  same 
monthly  ocean  thermal  climatology,  functioned  on  the 
same  schedule  ( one  analysis  per  day ).  and  had  access 
to  the  same  real-time  data  base  (i.e..  synoptic  ship, 
buoy,  bathy.  and  satellite  observations).  Throughout 
the  test.  OTIS  was  coupled  with  its  own  version  of 
TOPS  in  an  analysis-prediction-anaiysis  data  assim¬ 
ilation  cycle  independent  of  TEOTS.  The  TOPS  mixed- 
layer  model  coupled  with  OTIS  contained  exactly  the 
same  physics  and  used  exactly  the  same  atmospheric 
forcing  as  that  which  was  coupled  to  TEOTS. 

The  purpose  of  this  paper  is  to  document  the  tech¬ 
niques  and  assumptions  used  in  OTIS,  and  present  the 
results  of  the  operational  test  of  the  global-scale  OTIS 
at  FNOC. 


In  performing  an  analysis.  OTIS  proceeds  point-by¬ 
point  through  its  horizontal  grid,  producing  a  complete 
profile  from  the  surface  to  the  deepest  analyzed  level 
at  each  gridpoint  before  advancing  to  the  next.  At  each 
of  these  gridpoiats.  the  analy  sis  sequence  is  1  )  sea  sur¬ 
face  temperature.  2)  temperature  in  the  mixed  layer, 
and  3)  temperature  below  the  mixed  layer. 


OTIS  represents  the  analyzed  sea  surface  tempera¬ 
ture  ( SS  I  )  at  the  A  til  gridpoint  '/'*■'  as 


2.  Description  of  the  analysis  procedure 


a  Sea  surface  temperature  analysts 
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these  data  and  tend  to  control  their  weighting  relative 
to  climatology,  other  data,  and  TOPS  in  the  analysis. 
The  MCSST  data  arc  handled  in  a  special  way  which 
complicates  the  determination  of  (ff/’J^cssr.  Before 
assimilation  into  the  analysis,  the  MCSST  reports 
within  one-half  mesh  length  of  each  gridpoint  arc  block 
averaged  around  the  gridpoint  to  form  “super-obser¬ 
vations"  (see  DiMcgo  1 988)  assumed  valid  at  the  grid- 
point.  Because  high  MCSST  data  densities  usually  re¬ 
sult  in  many  redundant  MCSST  reports  around  a  grid- 
point,  this  averaging  process  provides  an  efficient  way 
to  incorporate  these  high-resolution  data  without  losing 
any  information  on  the  resolvable  thermal  field.  Thus, 
if  the  /  th  observation  in  Eq.  ( 1 )  is  an  MCSST  obser¬ 
vation,  it  is  in  fact  given  by 


nn 


MCSST 


A/ 


M 

V 

t-J 

y- 1 


(8) 


where  t,  is  the y'th  individual  MCSST  report  and  M  is 
the  total  number  of  individual  MCSST  reports  used  in 
the  block  average.  Thus,  in  view  of  Eq.  (8),  the  expres¬ 
sion  for  ( <r,")  mcsst  Becomes 


(Ok 


MCSST  ~ 


M 


(  <  klCSS  I  +  (O2  +  (O*]  +  C  MCSST 


(9) 


where  tMcssr  and  Cmcsst  arc  the  standard  deviation 
and  mean,  respectively,  of  the  instrumental  error  for 
the  MCSST  reports  and  o'  is  the  standard  deviation 
ofthcSST  field  resulting  from  resolvable  SST  variations 
over  the  domain  in  which  the  MCSST  block  averaging 
is  done.  The  term  ( a,')2  represents  the  variance  added 
because  the  differencing  of  the  MCSST  data  with  cli¬ 
matology  is  done  after  the  block  averaging  and  with 
the  climatological  temperature  valid  at  the  gridpoint 
rather  than  before  the  block  averaging  and  with  tem¬ 
peratures  obtained  by  interpolating  climatology  to  each 
individual  MCSST  observation.  This  term  is  estimated 
by  assuming  a  constant  reference  SST  gradient  and 
averaging  the  resulting  temperature  variance  over  one 
grid  space.  Thus, 


(o,'\ 3 


(TAX)1 

12 


(10) 


where  T’  is  a  reference  climatological  SST  gradient 
(taken  to  be  0.0 1°C  km-1 )  and  A, Vis  a  reference  grid 
space  for  the  OTIS  grid.  Thus,  for  the  ENOC  63  X  63 
hemispheric  grid  (AX  =  320  km),  o[  ~  ().85°C, 

No  attempt  is  made  to  correct  for  any  mean  error 
(i.c..  bias)  in  the  MCSST  data.  Rather,  (’mcsst  effec¬ 
tively  becomes  a  lower  bound  on  the  MCSST  obser¬ 
vational  error  when  M  becomes  large.  The  quantities 
(mc  sst  and  < 'mcsst  arc  user-specified  parameters  which 
lend  to  control  the  relative  weighting  of  the  MCSST 
data  in  the  analysis. 


The  value  of  the  TOPS  prediction  rms  error  of  for 
use  in  Eq.  (5)  is  obtained  from 

(  ffF)~  =  (ff*")  INITIAL  +  6/V  (II) 

where  Gu"  Initial  is  the  rms  error  of  the  resolvable 
thermal  field  from  the  previous  analysis  initializing 
TOPS  and  Gk  is  the  rms  TOPS  error  growth  over  one 
analysis -prediction-analysis  cycle  (i.c.,  24  hours). 
Thus,  following  the  approach  of  Benglsson  and  Gus- 
tafsson  i  1972),  the  error  in  the  predicted  thermal  field 
is  taken  to  be  the  error  in  the  initial  conditions  plus 
the  error  added  due  to  error  growth  in  the  prediction 
model.  The  quantity  G  is  a  user-specified  parameter 
chosen  to  reflect  the  rate  of  error  growth  in  TOPS.  It 
tends  to  control  the  weight  assigned  to  the  TOPS  pre¬ 
diction  relative  to  climatology  and  observations. 

Clancy  et  al.  ( 1989)  present  idealized  data  assimi¬ 
lation  studies  to  illustrate  the  behavior  of  the  analysis- 
forecast- anab  .is  coupling  between  OTIS  and  TOPS 
described  above.  In  general,  the  mixed-layer  model 
predictions  reduce  the  error  level  of  the  analysis  by 
accounting  for  the  response  of  the  upper  ocean  to  local 
atmospheric  forcing  and  carrying  the  integrated  effect 
of  this  response  forward  in  lime.  In  the  absence  of  large 
TOPS-predictcd  changes  in  the  thermal  structure  and 
observational  data  to  sustain  OTIS-analyzcd  thermal 
anomalies,  the  thermal  fields  decay  toward  climatology 
at  a  rate  which  is  governed  by  the  prescribed  TOPS 
error  growth  rate  Gk . 

The  last  user-specified  parameter  is  the  noise-to-sig- 
nal  ratio  for  the  subgrid-scale  error  X/'  defined  by 


Like  A.  B ,  and  C.  this  quantity  must  be  chosen  to 
reflect  the  statistics  of  ocean  thermal  variability.  For 
example.  X/'  must  be  larger  in  dynamically  active  re¬ 
gions,  where  subgrid-scalc  eddies  produce  large-am¬ 
plitude  but  unresolved  perturbations  in  the  thermal 
field,  than  in  more  quiescent  regions  where  the  eddy 
field  is  weak.  Of  course,  X,"  is  a  function  of  the  grid 
mesh  length,  becoming  smaller  as  the  grid  becomes 
finer  and  the  unresolved  portion  of  the  thermal  field 
decreases. 

Equation  ( 12)  is  used  along  with  the  rms  difference 
between  balhy  observed  temperatures  and  climatology 
<t/’  to  obtain  of  and  a/1.  Neglecting  (rathy-  which  is 
small  compared  to  o,”  for  the  grid  currently  utilized  by 
OTIS,  and  assuming  that  rff  and  <r”  arc  uncorrelated, 
(a,1’)2  can  be  written 

(«T,Y  =  (T,f)J  +  W)J  (13) 


or  just  the  sum  of  the  mean-square  resolvable  signal 
about  climatology  plus  the  mean-square  subgrid-scale 
noise.  From  (12)  and  (13), 


(<r,r)2 


I  +  X/’ 


(14) 
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(ff(")’  =  — - (  rr/’)-.  (15) 

'  1  +  X,"  ' 

It  is  not  possible  to  calculate  (a,1')1  directly  from  the 
available  real-time  bathy  data  at  every  observational 
location  simply  because  the  bathy  data  coverage  and 
observational  density  will  not  permit  it.  Therefore. 

( <x/’):  is  estimated  from 

,  14  X/'  ,,  , 

(  of  V  =  -  ( n,n U  (  16) 

*  '  'Nft 

where  ( rrfj, ):  is  the  mean-square  departure  of  the  bathy 
observations  from  climatology  and  X",  is  the  average 
value  of  X,",  both  calculated  over  the  entire  analysis 
domain.  The  data  used  in  this  calculation  consist  of 
a  i uniting  60-day  window  of  bathy  observations  re¬ 
ported  to  FNOC  in  real-time,  amounting  to  about  0000 
reports  globally. 

Using  the  assumption  of  Eq.  (16).  Eqs.  (14)  and 
(15)  become 


\  " 

(«/•')-  ^~r:  ):-  (IS) 

I  +  X" 

With  X,"  specified  by  the  user,  and  and  X",  calcu¬ 
lated  internally  b>  OIlS.  Eqs.  (17)  and  (IS)  provide 
n,  and  a,"  for  use  in  Eqs.  ( 4 ).  (  5 ).  ( 6 ).  (  7 ).  ( 0 ).  and 
(12).  Note  that  .  and  thus  n;  and  of. are  determined 
separately  for  each  standard  level  of  the  OHS  vertical 
output  grid  (see  section  2b).  Note  also  from  Eq.  (17) 
that  it,'  is  independent  of  horizontal  location.  Thus, 
implicit  in  Fq.  (  16)  is  the  assumption  that  variations 
of  the  rms  bathy  error  nf  with  location  arc  produced 
entirely  by  variations  in  the  rms  subgrid-scale  error  n 
Equations  ( 2 )— ( 12)  and  { 1 7 )— (  18)  constitute  a 
closed  set  which  can  be  solved  for  the  weights  <\k,  and 
dt .  Once  these  weights  are  obtained,  the  analyzed  tem¬ 
perature  is  calculated  from  Eq.  (  I  ).  Also,  the  rms  error 
of  the  resolvable  thermal  field  nk"  can  then  be  calculated 
from 

\ 

(itUU  <0:(l  fa  -  (19) 

i  i 

(see  White  and  Bernstein  1474).  Recognizing  that 
OTIS  initializes  the  TOPS  forecast  which  will  provide 
//'  for  the  following  day's  analysis,  the  value  of  (  a k“): 
obtained  from  (  19)  becomes  ( 'TA“)i\m\i  of  equation 
(  I  I  )  ( i.e..  the  error  in  the  TOPS  initial  conditions )  for 
the  following  day’s  analysis. 

Finally,  the  expected  departure  of  the  thermal  field 
from  an  independent  bathy  observation  made  precisely 
at  the  gridpoint  at  analysis  time  ok  is  given  by 


with  the  three  terms  inside  the  brackets  representing 
the  contributions  of  analysis  error,  subgrid-scale  error, 
and  instrumental  error. 

h  Mixed-layer  analysis 

Upon  conclusion  of  the  SST  analysis  at  gridpoint  A. 
OHS  performs  an  Ol  analysis  for  mixed-layer  depth 
(  MED).  In  analogy  to  Eq.  (  I  ).  the  expression  for  the 
analyzed  MED  hk“  is 

V 

lh“  =  hk  +  <u,(/f,"  -  li,' )  4  fikUhr  -  lik  )  (21  ) 

;  I 

where  hk  is  the  climatological  MED  at  the  gridpoint. 
/;/’  the  predicted  MED  at  the  gridpoint  from  the  24- 
hour  TOPS  forecast  from  the  previous  day’s  analysis. 
/;/'  the  observed  MED  at  location  /.  h,'  the  climatolog¬ 
ical  MED  at  location  ;.  <» y,  the  weight  applied  to  the 
/ th  observed  MED  anomaly  (//,"  -  //,' )  assimilated  at 
the  gridpoint.  dk  the  weight  applied  to  the  I  OPS  pre¬ 
dicted  MED  anomaly  (  hkr  -  hk )  at  the  gridpoint.  and 
V  the  number  of  observations  assimilated  at  the  grid- 
point. 

The  climatological  MED//f  is  determined  from  ex¬ 
amination  of  the  climatological  temperature  profile, 
defined  at  gridpoint  A  on  the  OTIS  vertical  grid  (see 
Table  I  ).  Beginning  at  the  surface  and  working  down¬ 
ward.  vertical  temperature  gradients  are  examined  to 
find  the  shallowest  pair  of  levels  between  which  the 
vertical  gradient  is  less  than  0.05 °C '  m  '.  The  layer 
depth  hk  is  then  taken  simply  as  the  depth  of  the  shal¬ 
lower  level  of  this  pair,  and  is  thus  alway  s  one  of  the 
discrete  levels  listed  in  Table  1.  Exactly  the  same  pro¬ 
cess  is  carried  out  on  the  TOPS  and  bathy  profiles  to 
define  hkr  and  h,".  The  TOPS-predicted  profile  is  de¬ 
fined  on  the  same  vertical  grid  as  OTIS,  but  the  depths 
defining  the  bathy  profiles  are  arbitrary,  and  thus  the 
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h,'” s  generally  fall  somewhere  in  between  the  OTIS/ 
TOPS  fixed  levels  of  Table  1. 

The  weights  «*,  and  /J*  in  Eq.  (  21 )  are  obtained  by 
solving  the  N  +  1  equations  of  ( 2 )  as  before.  In  solving 
(2),  the  parameters  ?/„,  X,"  and  X/’  are  assumed  pro¬ 
portional  to  the  corresponding  values  for  the  temper¬ 
ature  analysis  at  100  m  depth.  This  implies 

(<b°)toi.D  =  ( A )(  (t/')100  m  (22) 

(<7*P)mlD  =  (A)(ff7)lOO  m  (23) 

(<Ta‘)mlD  =  (X)(ff(‘)ioo  m  (24) 

where  the  subscript  MLD  indicates  values  are  for  MLD 
variations,  and  K  is  a  constant. 

Although  hk  and  hkp  can  only  take  on  the  discrete 
values  of  Table  I,  hku  resulting  from  the  solution  of 
(21 )  and  (2)  can  take  on  any  value.  Thus,  it  can  be 
thought  of  as  a  “floating  level"  in  the  analysis  system. 

With  the  analysis  for  hk“  complete,  OTIS  sets  the 
temperatures  at  the  gridpoints  of  Table  1  which  are 
shallower  than  hk“  according  to 

»] 

(  1\“  )m  =  (  7'a‘,)<,  +  I  ( A '//’)„  (25) 

a  i 

where  (  7*." )„,  is  the  analyzed  temperature  at  level 
(  7\“ )0  is  the  analyzed  SST.  and  (  A//),,  is  the  change 
in  TOPS  predicted  temperature  from  level  n  -  I  to 
level  n  (level  0  implies  the  surface).  Thus,  the  shape 
of  the  temperature  profile  from  the  surface  to  the  base 
of  the  mixed  layer  (e.g..  isothermal,  weakly  stratified, 
or  multiple  thermocline )  is  controlled  exactly  by  TOPS. 

< :  Siib-mi.xeii-layer  analysis 

Following  conclusion  of  the  mixed-layer  analysis  at 
gridpoint  A;.  OTIS  performs  an  analysis  for  temperature 
at  and  below  the  analyzed  MI  D.  Below  this  level.  OTIS 
utilizes  a  completely  variable  or  “floating"  vertical  grid 
designed  to  concentrate  resolution  in  regions  of  high 
vertical  gradient.  Beginning  at  the  MLD  and  working 
downward,  the  mesh  spacing  of  this  floating  grid  Ac  is 
defined  as 

Ac  50  m  t  (400  m:C  1  )•-  — 

<): 


for  5  m  §  A:  ^  SO  m  (  26 ) 

with  r  taken  positive  downward  from  the  sea  surface. 

With  the  vertical  grid  extending  from  the  base  of 
the  mixed  layer  downward  to  400  m  defined  by  (  26  ), 
OTIS  performs  an  analysis  for  temperature  at  each  of 
these  levels.  This  analysis  is  carried  out  by  solving  (  1  ) 
and  ( 2 )  at  each  of  the  floating  levels. 


In  general,  the  bathy  observations  which  contribute 
to  the  sub-mixed-layer  analysis  extend  to  varying 
depths.  Bathys  which  reach  to  less  than  250  m,  or  which 

have  —  1(7/’-  7/  )  |  >0.01°C  m  1  near  their  deepest 
a: 

reported  depth  zB,  are  used  in  the  analysis  only  for 
levels  above  zB.  Bathys  extending  to  greater  than  250 
m  depth  for  which 

A 

—  1(77  -  77)1  <  o.orc  m'1 

oz 

near  zB  are  used  in  the  analysis  at  all  levels  down  to 
400  m.  however,  by  extrapolating  ( T,"  -  T‘  )  down¬ 
ward.  The  mean-square  error  added  to  the  extrapolated 
anomaly  below  zB  by  this  process  is  estimated  from 

(<r/*)f.x  =  10.0225  +  0.02 

X  exp[  -0.05 ( r*  -  250)] } ( r  -  z„)  ( 27 ) 

where  r  is  the  depth  of  the  extrapolation  point  for  the 
ith  bathy  in  meters.  Equation  (7)  for  extrapolated 
bathy  observations  below  zH  then  becomes 

(<T/’)fWIIIY  =  <BMMY  +  (ff/')‘  +  (<I|")fcX-  (28) 

The  quantities  i)„,  X/’.  and  X/'  required  in  (2)  are 
derived  from  the  user-specified  parameters  .  L  .  Bk .  C\. 
X/'  and  «b.mhy  as  before.  In  calculating  X,"  arid  X/’  via 
(4).  (5).  and  (7).  the  quantities  a,',  a,"  and  a,r  are 
interpolated  linearly  from  the  fixed-levels  of  Table  1 
to  the  floating  levels.  Equation  ( 28  )  is  used  in  place  of 
( 7  )  for  extrapolated  bathy  anomalies  below  r/(  as  dis¬ 
cussed  above.  Note  that  the  inverse  decorrelation  scales 
b.  Hk .  C\.  and  the  noise-to-signal  ratio  X,"  for  sub- 
mixed-layer  analysis  can  be  assigned  v  alues  by  the  user 
distinctly  different  from  those  at  the  surface.  This  is 
particularly  important  for  C .  which  generally  has  a 
much  smaller  value  below  the  MLD  reflecting  the 
much  longer  time  scale  variability  there. 

Upon  conclusion  of  the  sub-mixed-layer  tempera¬ 
ture  analysis,  the  analyzed  temperature  anomalies  at 
the  floating  levels  defined  by  (26)  are  vertically  inter¬ 
polated  to  the  fixed  levels  of  Table  I  between  the  MLD 
and  400  m  depth,  added  back  to  climatology  ,  and 
output. 

3.  User-specified  parameters 

The  nine  user-specified  parameters  used  in  O  US  are 
listed  in  fable  2.  The  values  of  <sm,.  and  (iiwm  are 
obtained  from  Earle  (  1085)  and  White  and  Bernstein 
(  1070  ).  respectively  .  The  values  of  <M<  SSi  and  <\u  sm 
arc  based  on  Hawkins  cl  al.  (  1086)  and  Strong  and 
McClain  (  1085).  I  he  values  of  I  and  li  range  from 
1/(2000  km)  in  low-latitude  open-ocean  regions  to 
I  /( 200  km  )  in  coastal  regions.  The  values  oft/  and  X" 
also  vary  widely .  In  general,  (i  is  smaller  in  midlatitude 
open-ocean  regions  than  in  tropical  or  western  bound- 
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Table  2.  User-specified  parameters. 


Parameter 

Value 

•SHIP 

2.29°C 

•bathv 

0.20°C 

•mcsst 

0.7  l°C 

^MCSST 

0.22-C 

c 

1/(20  days)  above  the  MLD.  1/(60  days) 
below  the  MLD. 

.1.  B.  O. 

Variable,  prescribed  al  the  center  of  10°  x  20° 
latitude/longitude  rectangles  and  interpolated 
linearly  to  gridpoint  locations.  Separate  values 
for  above  and  below  the  MLD. .  t  ranges  from 
1/(200  km)  to  1/(2000  km),  B  ranges  from  1/ 
(200  km)  to  l/(  1000  km),  (i  ranges  from  0. 1 0 
to  l.0°C,  and  A"  ranges  from  0.8  to  6.0. 

ary  current  regions,  where  the  atmospheric  forcing  and / 
or  physics  of  TOPS  is  not  fully  adequate.  In  general, 
X"  is  smaller  in  quiescent  regions  ( e.g.,  eastern  Pacific ) 
than  in  dynamically  active  regions  characterized  by 
strong  mesoscale  eddies  which  are  unresolved  by  the 
horizontal  grid  (e.g..  Kuroshio  or  Gulf  Stream  regions). 

4.  Data  selection  and  screening 

Before  the  analysis  begins,  thinning  and  quality- 
control  algorithms  are  applied  to  the  data.  First,  the 
bathy  data  are  screened  by  applying  a  sorting  process 
which  limits  the  bathy  data  density  in  2.5°  X  2.5°  lat¬ 
itude/longitude  squares  to  the  60  most  recent  reports, 
with  no  more  than  30  taken  from  the  same  platform. 
This  thinning  process  is  a  practical  necessity  for  han¬ 
dling  ocean  weather  stations,  which  typically  report 
two  bathy  observations  per  day. 

Prior  to  computing  the  climatological  anomalies 
from  the  observations,  all  data  are  checked  for  physi¬ 
cally  unreasonable  temperatures.  Any  observations  re¬ 
porting  less  than  -2°C  or  more  than  40°C  are  dis¬ 
carded.  Once  the  observed  anomalies  are  computed 
and  passed  to  the  analysis,  they  are  subjected  to  a  gross 
error  check.  Any  observations  which  deviate  by  more 
than  7°C  from  climatology  are  rejected. 

In  addition,  OTIS  applies  a  ship-tracking  algorithm 
to  screen  out  bathy  reports  with  obvious  position  errors. 
These  types  of  errors  are  fairly  common  in  the  real¬ 
time  bathy  dataset  and  result  primarily  from  coding 
and  radio  transmission  errors  (e.g.,  latitude  and  lon¬ 
gitude  transposed,  wrong  hemisphere,  etc.).  Thus, 
OTIS  monitors  the  successive  positions  of  each  ship 
reporting  bathys  at  any  time  during  the  previous  60- 
dav  period  and  flags  jumps  in  position  which  would 
require  impossibly  large  ship  speeds  (bathys  deployed 
from  aircraft  are  not  tested  in  this  manner).  Bathy 
observations  corresponding  to  the  erroneous  position 
reports  are  excluded  from  further  processing  by  OTIS. 

For  the  analysis  at  each  gridpoint.  OTIS  utilizes  only 
data  which  fall  within  certain  temporal  and  spatial 
windows  (see  Phoebus  1988).  For  the  SST  analysis. 


the  temporal  window  extends  from  the  OTIS  analysis 
time  back  60  hours.  For  the  MLD  and  sub-mixed- 
layer  analyses,  the  window  extends  back  60  days.  The 
spatial  window  is  defined  in  terms  of  the  spatial  cor¬ 
relation  scales  I j  A  and  \/B.  The  window  is  defined 
as  plus  or  minus  one  correlation  scale  in  the  north- 
south  and  east-west  directions  from  the  gridpoint.  At 
gridpoints  where  the  east-west  correlation  scale  is  larger 
than  the  north-south  correlation  scale,  as  is  often  the 
case,  this  results  in  a  spatial  window  that  is  rectangular 
in  shape,  extending  farther  in  the  east-west  direction 
than  in  the  north-south  direction.  Variations  of  A  and 
B  w  ith  location  will  alter  the  size,  shape,  and  orientation 
of  this  window  from  gridpoint  to  gridpoint.  and  at  var¬ 
ious  depths  for  the  same  grid  location. 

The  area  within  the  appropriate  spatial  window  is 
searched  and,  ideally  ,  all  observations  within  the  spec¬ 
ified  window  are  collected.  In  practice,  however,  the 
number  of  observations  that  can  be  collected  is  limited 
to  700  ( a  limit  which  is  seldom  exceeded ).  During  the 
data  search,  the  correlation  of  each  observation  with 
the  gridpoint  r),k  is  computed  from  Eq.  (3).  Once  all 
of  the  observations  are  collected,  they  are  sorted  by 
these  correlations  to  select  the  1 5  reports  most  highly 
correlated  w  ith  the  gridpoint  location  for  potential  use 
in  the  analysis.  If  fewer  than  1 5  observations  are  found, 
all  are  retained,  subject  to  one  final  quality  control 
procedure. 

The  ( up  to )  15  selected  reports  are  screened  for  spu¬ 
rious  data  using  a  horizontal  consistency  check,  or 
buddy  check,  patterned  after  that  of  DiMego  ct  al. 

( 1985 ).  The  particular  details  of  the  OTIS  buddy  check 
are  described  by  Phoebus  ( 1 989 ).  Briefly  ,  the  inequality 

| (  T;,  _  yy  )  _  (  yy  _  yy  )|  (a  -  hVll ) <rk'  ( 29 ) 

is  examined  for  each  pair  of  observations.  If  the  in¬ 
equality  is  satisfied,  observations  /  and  j  corroborate 
each  other  and  arc  thus  both  given  “keep  flags"  equal 
to  the  correlation  between  them  9,,:  if  the  inequality  is 
not  satisfied,  the  observations  contradict  each  other 
and  are  thus  both  given  “toss  flags"  equal  to  9,,.  Fol¬ 
lowing  examination  of  all  remaining  pairs  of  reports, 
the  single  observation  with  the  highest  summation  of 
toss  flags  is  removed  if  the  summation  exceeds  2.  How¬ 
ever.  observations  whose  keep  flags  sum  to  2  or  more 
are  retained  regardless  of  toss  flag  settings.  This  process, 
beginning  with  the  examination  of  the  inequality  ( 29) 
and  with  all  flag  summations  reset  to  0.  is  then  repeated 
until  no  observations  have  toss  flags  which  sum  to  more 
than  2.  The  remaining  observations  are  utilized  in  the 
analysis  at  the  gridpoint. 

Equation  (29)  states  that  the  absolute  temperature 
ditference  between  two  anomalies  is  compared  to  some 
tolerance  which  is  a  function  of  the  correlation  between 
them  i)„  and  the  long-term  rms  variation  of  the  thermal 
anomalies  about  climatology  <t<'.  Conceptually,  obser¬ 
vations  which  are  more  highly  correlated  with  one  an¬ 
other  are  expected  to  agree  more  closely,  and  temper- 
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ature  differences  in  regions  of  high  expected  rms  de¬ 
partures  from  climatology  are  given  more  tolerance. 
The  use  of  the  correlation  »?„  as  the  flag  value  is  also  a 
subtle  way  to  control  the  impact  one  observation  has 
on  the  retention  or  rejection  of  another.  Observations 
which  are  highly  correlated  contribute  more  to  the  toss 
or  keep  Hag  summations  of  one  another  than  do  those 
observations  which  are  farther  separated  in  space  and 
time.  Currently,  the  values  of  a  and  b  are  set  to  3.0 
and  1 .5,  respectively.  These  values  were  chosen  simply 
through  trial  and  error  (see  Phoebus  1980). 

5.  Qualitative  comparison  of  OTIS  and  TKOTS 

a.  Surface  thermal  .fields 

The  OTIS  SST  field  for  28  March  1988  is  show  n  in 
Fig.  1 .  This  graphic,  and  the  ones  to  follow,  were  pro¬ 
duced  by  interpolating  the  fields  (via  a  Bessel  tech¬ 
nique)  from  the  standard  FNOC  NHEM  and  SHEM 
63  X  63  polar  stereographic  grids  to  the  standard  FNOC 
73  X  144.  2.5  degree  spherical  grid,  and  then  con¬ 
touring. 

Though  the  grid  resolution  in  the  model  is  too  coarse 
to  accurately  represent  fronts,  OTIS  shows  large  hor¬ 
izontal  SST  gradients  in  the  South  China  Sea.  the  Sea 
of  Japan,  the  Kuroshio  region  and.  particularly,  in  the 
Gulf  Stream  region.  The  tropical  western  Pacific  SST 
field  is  quite  featureless,  w  hile  the  eastern  tropical  Pa¬ 
cific  shows  much  structure  suggestive  of  upwelling 
along  the  equator  and  off  the  coast  of  Peru. 

The  OTIS  SST  field  for  28  April  1988  is  shown  in 
Fig.  2.  Comparison  of  Figs.  1  and  2  illustrates  the  ex¬ 
pected  temporal  variability  of  the  large-scale  thermal 
field  over  one  month's  time;  OTIS  shows  SST  increases 


in  the  upper-  and  midlatitudes  of  the  Northern  Flemi- 
sphere  associated  with  the  springtime  warming.  A 
much  more  prominent  change  over  this  period,  how¬ 
ever.  is  the  cooling  along  the  equator  in  the  eastern 
Pacific.  Most  likely  due  to  enhanced  equatorial  up- 
welling  in  response  to  stronger  than  normal  easterly 
w  inds  associated  with  the  "La  Nina"  weather  pattern, 
this  cooling  trend  led  to  the  coldest  June  SST  anomalies 
in  the  eastern  tropical  Pacific  since  the  early  1970s 
(Climate  Analvsis  Center  1988). 

The  OTIS  SST  anomaly  field  for  28  March  1988  is 
shown  in  Fig.  3.  Numerous  features  are  prominent, 
including;  1  )  a  cold  anomaly  extending  across  most  of 
the  Pacific  between  35°  and  45° N.  2)a  warm  anomaly 
oriented  WSW-ENE  in  the  subtropical  western  Pacific. 
3)  a  warm  anomaly  in  the  eastern  North  Pacific  ex¬ 
tending  southwestward  from  Baja  California.  4 )  a  cold 
anomaly  along  the  equator  at  1 25°W.  5 )  a  cold  anom¬ 
aly  in  the  South  Pacific  centered  at  40°S.  1 10° W.  6) 
a  warm  anomaly  off  the  coast  of  Chile  at  30°S,  90°  W. 
and  7)  a  warm  anomah  in  the  South  Atlantic  at 
2()°S.  1 7°W. 

The  OTIS  SST  anomaly  field  on  28  April  1988  is 
show  n  in  Fig.  4.  Comparison  of  Figs.  3  and  4  illustrates 
changes  in  the  OTIS  SST  anomaly  field  over  one 
month's  time.  The  midlatitude  cold  anomaly  across 
most  of  the  Pacific  on  28  March  extends  farther  south 
and  shows  a  larger  amplitude  on  28  April.  The  WSW- 
ENE  warm  anomaly  in  the  subtropical  western  Pacific 
on  28  March  is  essentially  gone  bv  28  April,  replaced 
by  two  small  warm  anomalies  centered  at  1 5°N.  1 32°E 
and  22°N.  172°W.  and  a  larger  warm  anomaly  cen¬ 
tered  at  26 °N.  I60°E.  The  warm  anomaly  extending 
southwestward  from  Baja  California  on  28  March  ex- 


f-'Ui .  I .  Sea  surface  temperature  (SSI)  field  from  OTIS  on  78  March  1 088.  The  contour  interval  is  I  ('. 
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Lie;.  2.  Sea  surface  temperature  ( SS  I  )  field  from  O I  IS  on  28  April  1488  I'he  contour  interval  is  I 


tends  farther  north  and  west  by  28  April,  and  achieves 
a  secondary  maximum  ENE  of  Hawaii.  The  cold 
anomaly  along  the  equator  at  125°W  on  28  March 
moves  westward  to  about  1 35° W  by  28  April,  and  an 
even  more  intense  equatorial  cold  anomaly  forms  be¬ 
tween  85°  and  1()5°W.  The  cold  anomaly  at  40°S. 
1 10° W  on  28  March  weakens  and  moves  northwest¬ 
ward  to  35°S.  ;20°W  bv  28  April.  The  warm  anomalv 
off  the  coast  of  Chile  at  3()°S.  90°  W  on  28  March  re¬ 


mains  evident  but  weakened  on  28  April.  The  warm 
anomaly  in  the  South  Atlantic  at  3()°S.  17°W  on  28 
March  weakens  and  moves  westward  to  30°S.  27°\V 
by  28  April. 

The  TEOTS  SST  and  SST  anomaly  fields  on  28 
March  and  28  April  (not  shown)  exhibit  very  little 
qualitative  difference  with  the  OTIS  fields  of  Figs.  1- 
4.  Thus,  though  they  employ  rather  different  data  as¬ 
similation  techniques  (Hoi!  et  al.  1979;  Clancy  and 


Ik.  L  Sea  surface  temperature  (  SS I  )  anomalv  field  from  O' t  IS  on  28  March  1488  I  he  contour 'interval  is  U  S  ( 
Positive  contours  are  dashed  and  negative  contours  are  displaced  with  a  thin  solid  line 
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I  Hi.  4  Sea  surface  temperature  (SSI)  anomalv  field  from  OHS  on  28  April  ll*SS.  The  contour  inters  al  is  O.V( 
Positive  contours  are  dashed  and  negative  contours  arc  displaced  with  a  thin  solid  line. 


Poliak  1983).  the  two  analysis  systems  produce  very 
similar  SST  products  on  the  global  scale.  This  probably 
retlec's  the  abundance  of  MCSST  data  available  to 
them. 

Also,  the  28  March  and  28  April  OT  IS  SST  anomaly 
fields  arc  generally  in  good  agreement  with  the  SSI 
anomaly  fields  of  NO  A  A  (  1988a.  1988b.  1988c).  Note 
that  a  direct  quantitative  comparison  of  the  OT  IS  and 
NO.AA  SST  anomaly  fields  is  beyond  the  scope  of  this 
paper  and  is  not  particularly  meaningful  in  any  case 


since  the  NOAA  product  represents  a  monthly  aver¬ 
aged  field  while  the  OT  IS  product  represents  a  daily- 
averaged  field. 

/>.  Suh.surf  uc  thermal  hclc/s 

The  OT  IS  and  PLOTS  temperature  al  200  m  depth 
(T200)  fields  for  28  March  1988  arc  shown  in  Fig.  5. 
In  general.  OHS  produces  sharper  horizontal  gradients 
in  the  1200  thermal  field,  particularly  in  the  subtropical 


I  k.  5a.  Temperature  al  200  m  depth  (1200)  from  OT  IS  on  28  March  1488.  The  contour  mterv...  is  1C. 
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I  k;.  5b.  As  in  l  ig  5a.  bui  lor  I  LOTS. 


Pacific  and  Atlantic.  OHS  provides  a  better  represen¬ 
tation  of  the  western  boundary  currents  at  200  m  depth 
than  I  HOTS,  including  I  )  separation  of  the  high  gra¬ 
dient  /ones  of  the  Kuroshio  and  the  Oyashio.  2  )  a  more 
physically  realistic  orientation  and  packing  of  the  iso¬ 
therms  in  the  Kuroshio  south  of  Japan.  3 )  tighter 
packing  of  isotherms  in  the  Gulf  Stream  extension  re¬ 
gion  between  60°  and  35°W.  and  4  )  a  more  physically 
realistic  orientation  and  packing  of  isotherms  west  of 
70°W  and  south  of  Cape  Hatteras.  In  addition,  the 


OTIS  T200  field  shows  a  clear  ( though  poorly  resolved ) 
representation  of  the  Loop  Current  in  the  Gulf  of 
Mexico.  This  feature  is  all  hut  absent  in  the  PHOTS 
T200  field. 

The  OTIS  and  TH.OTS  temperature  anomaly  fields 
at  200  m  depth  (  T  200  anomaly)  for  2K  March  llJ88 
are  shown  in  Fig.  6.  Substantial  differences  are  apparent 
in  the  T  200  anomaly  fields  produced  by  the  two  mod¬ 
els.  In  General,  the  PHOTS  T200  anomaly  field  appears 
noisii  nan  the  corresponding  OTIS  field,  particularly 


I  k,  6a.  I empcraturc  anomaly  al  200  m  depth  i  I  200  anomalv  )  from  Ol  IS  mi  2S  March  loss  I  he  contour 
interval  is  0.5  C  Positive  contours  arc  dashed  anil  negative  contours  arc  display cvl  with  a  thin  solid  line 
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in  the  tropical  Pacific.  The  features  in  the  OTIS  anom¬ 
aly  field  generally  show  more  spatial  continuity  and 
exhibit  a  more  physically  reasonable  east-west  (vice 
north-south  or  circular)  pattern.  A  large  warm  anom¬ 
aly  extending  across  most  of  the  tropical  South  Pacific 
between  5°  and  25°S  is  present  in  both  the  O  VIS  and 
TEOTS  fields,  but  the  northern  and  southern  bound¬ 
aries  of  this  anomaly  are  represented  much  more 
sharply  by  OTIS.  This  corresponds  to  the  sharper  hor¬ 
izontal  gradients  in  the  OT  IS  T200  field  for  this  region 
which  were  obvious  in  Fig.  5. 

In  some  regions,  the  differences  between  the  OTIS 
and  TEOTS  1200  anomalies  are  quite  large.  For  ex¬ 
ample.  TEOTS  shows  a  w  arm  anomaly  of  over  2.0°C 
at  47°N.  177°W.  but  this  feature  is  not  present  in  the 
OTIS  anomaly  field.  An  even  more  striking  case  is  the 
TEOTS  warm  anomaly  of  over  4.0°C  olf  the  coast  of 
New  Foundland  at  45°N.  4>°W:  OTIS  shows  a  cold 
anomaly  with  amplitude  greater  than  l.5°C  in  this  re¬ 
gion. 

I  he  large  differences  between  the  OTIS  and  TEOT  S 
T  200  anomaly  fields  arc  a  result  of  the  fact  that  TOTS/ 
TEOTS  performs  a  vertical  blending  process  from  the 
surface  to  400  m  depth  (  Moll  et  al.  1074:  C  lancy  and 
Poliak  19X3).  In  the  absence  of  adequate  bathy  data 
density,  this  tends  to  tie  subsurface  anomalies  (even 
below  the  mixed  layer)  to  the  SST  anomaly.  T  hus,  it 
is  never  apparent  whether  subsurface  anomalies  in 
EOTS/TEOTS  are  supported  by  bathy  data,  or  arc  the 
result  of  this  vertical  blending  process.  Also,  in  the  ab¬ 
sence  of  adequate  bathy  coverage.  EOTS/TEOTS  can 
produce  a  subsurface  thermal  field  which  is  worse  than 
climatology  if  the  real  thermal  anomaly  changes  sign 
with  depth,  l  or  example,  if  there  is  sufficient  surface 
data  to  map  a  warm  anomaly  at  the  surface,  but  in¬ 


sufficient  bathy  data  to  map  the  field  at  200  m  depth, 
then  the  EOTS  vertical  blending  will  tend  to  force  in 
a  warm  anomaly  at  200  m.  But.  if  there  is  actually  a 
cold  anomaly  at  200  m.  then  the  EOTS/TEOTS  tem¬ 
perature  at  this  depth  will  differ  from  reality  by  even 
more  than  climatology  (i.e..  it  will  be  warmer  than 
climatology  while  reality  is  colder  than  climatology  ). 
OTIS,  on  the  other  hand,  performs  its  equivalent  of 
the  EOTS/TEOTS  vertical  blending  process  only  from 
the  surface  down  to  the  base  of  the  mixed  layer.  Thus, 
temperature  anomalies  below  the  mixed  layer  in  OTIS 
can  arise  only  if  the  bathv  data  support  them. 

T  he  OTIS  and  TEOTS  T200  anomaly  fields  for  28 
April  1988  are  shown  in  Eig.  7.  As  was  the  case  for  28 
March,  substantial  differences  between  the  OTIS  and 
TEOTS  anomalies  are  evident.  OTIS  shows  warm 
anomalies  in  the  western  North  Atlantic,  roughly  along 
the  path  of  the  Gulf  Stream,  while  TEOTS  show  s  cold 
anomalies  in  this  region.  I'EOTS  shows  warm  anom¬ 
alies  in  the  Indian  Ocean  (  15°S.  9()°E),  the  eastern 
tropical  Pacific  (5°S-2()°N.  9()°- I40°VV ).  and  the 
western  tropical  Atlantic  (2()°S-I5°.N.  2()°-60°W) 
which  are  not  present  in  the  O  I  IS  field. 

Comparison  of  Figs.  6a  and  7a  illustrate  changes  in 
the  OT  IS  1200  anomaly  field  on  a  one-month  time 
scale.  Most  of  the  features  prominent  on  28  March  arc 
still  evident  on  28  April.  An  exception  is  the  cold 
anomaly  at  40°S.  I40°VV  on  28  March:  this  feature 
weakens  substantially  by  28  April.  Also,  the  1.5  degree 
warm  anomaly  southwest  of  Hawaii  (  17°N.  I65°W) 
on  28  March  increases  to  a  3.0  degree  anomaly  by  28 
April,  finally,  a  warm  anomaly  of  more  than  2  degrees 
forms  in  the  western  tropical  Pacific  at  4"N.  136  [ 
between  28  March  and  28  April. 

The  OT  IS  ami  TEOTS  temperature  anomaly  fields 
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Flo.  7a.  Temperature  anomaly  at  200  m  depth  (1200  anomaly  )  from  OTIS  on  2S  April  19SX.  The  eontour 
inter\al  is  0.5“C.  Positive  contours  are  dashed  and  negative  contours  are  displayed  with  a  thin  solid  line. 


at  400  m  depth  (T400  anomaly )  for  28  April  1988  are 
shown  in  Fig.  8.  Large  differences  in  these  fields,  which 
are  also  a  result  of  the  TFOTS  vertical  blending  dis¬ 
cussed  above,  are  apparent.  For  example,  in  response 
to  and  consistent  with  the  bathy  data  at  400  m  depth. 
OTIS  shows  a  cold  anomaly  in  the  central  South  Pacific 
between  10  and  20°S.  In  primary  response  to  its  own 
vertical  blending  process  and  without  the  support  of 
bathy  data  at  400  m  depth,  however.  TFOTS  shows  a 


warm  anomaly  in  this  region  and  throughout  most  of 
the  subtropical  South  Pacific. 

6.  Quantitative  verification  of  OTIS  and  TFOTS 

u.  Procedure 

Verification  statistics  for  OTIS  and  TEOTS  were 
compiled  from  comparison  against  independent  bathy 
data.  Although  the  bathy  data  are  subject  to  errors 


I- It  i .  7b.  As  in  F  ig.  7a.  hul  lor  I  LOTS. 
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Fig.  8a.  Temperature  anomaly  at  400  m  depth  (T400  anomaly)  from  OTIS  on  28  April  1488.  The  contour 
interval  is  0.5°C.  Positive  contours  are  dashed  and  negative  contours  are  displayed  with  a  thin  solid  line. 


which  can  be  quite  large  (including,  in  the  context  of 
model  comparison,  errors  produced  by  subgrid-scale 
features  unresolved  by  the  model  grids),  they  still  pro¬ 
vide  an  important  ground  truth  reference  for  ranking 
the  skill  of  ocean  thermal  models. 

Once  per  day  the  thermal  fields  from  OTIS  and 
TEOTS  are  interpolated  via  a  Bessel  technique  to  the 
locations  of  all  bathy  observations  received  at  FNOC 
during  the  previous  24  hours  but  not  yet  assimilated 


into  the  models  (making  the  observations  independent 
of  the  models).  Each  bathy  observation  is  then  verti¬ 
cally  interpolated  linearly  to  a  standard  vertical  grid 
between  the  surface  and  300  m.  and  spurious  obser¬ 
vations  are  discarded  automatically  by  a  bathy  error 
detection  algorithm,  which  is  less  sophisticated  than 
the  OTIS  quality  control  procedures  discussed  in  sec¬ 
tion  4. 

The  remaining  data  are  used  to  calculate  analyzed 


Fig.  8b.  As  in  Fig.  8a.  but  for  TEOTS. 
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minus  observed  temperatures  (i.e.,  the  "apparent  er¬ 
rors"  in  the  analyzed  temperatures)  for  both  OTIS  and 
TEOTS  at  the  standard  grid  depths  ofO.  12.5.  25.  50. 
75.  100.  150.  200.  and  300  m  at  each  baths  location. 
The  resulting  apparent  errors  for  both  OTIS  and 
IEOTS  from  each  individual  baths  at  each  standard 
depth  are  accumulated  in  a  running  90-das  history  file 
and  segregated  into  the  nine  verification  areas  shosvn 
in  Fig.  9.  From  the  data  in  this  histors  file,  the  apparent 
root-mean-square  (rms)  errors  for  both  OTIS  and 
TEOTS  are  calculated  at  the  standard  depths  in  each 
of  the  nine  verification  areas  from  the  90-day  accu¬ 
mulation  of  data  covering  the  months  of  February. 
March,  and  April  1988.  The  total  number  of  data  col¬ 
lected  and  used  in  these  calculations  amounts  to  9980 
bathv  observations. 

h  I  ertieal  profiles  of  apparent  error 

Plots  of  vertical  profiles  of  the  apparent  rms  errors 
for  OTIS  and  TEOTS  for  February -April  1988  for  each 
of  the  nine  verification  areas  are  shosvn  in  Figs.  lOa-i. 
The  shapes  of  the  error  profiles  for  OTIS  and  TEOTS 
in  each  area  are  quite  similar.  All  areas  except  the  mid¬ 
latitude  Atlantic  areas  (Areas  3  and  4.  Eigs.  10c  and 
lOd )  exhibit  local  maxima  in  the  error  profiles  at  some 
intermediate  depth  for  both  OTIS  and  TEOTS.  These 
local  maxima  are  quite  pronounced  for  the  Indian 
Ocean  (Area  5.  Fig.  lOe).  the  tropical  eastern  Pacific 
(Area  I  A.  Fig.  I  Of),  and  the  tropical  eastern  Atlantic 
(  Area  3A.  Fig.  lOh).  and  are  probably  due  primarily 
to  noise  in  the  verification  data  caused  by  unresolved 
vertical  displacements  of  the  strong  thermoelincs  in 
these  areas  produced  by  eddies  and  internal  waves.  The 
midlatitude  eastern  Atlantic  (  Area  3.  Fig.  I()c )  is  char¬ 
acterized  by  local  minima  in  the  error  profiles  for  both 


OTIS  and  TEOTS  at  about  100  m  depth.  The  midlat¬ 
itude  western  Atlantic  (Area  4.  Fig.  lOd)  shows  a 
monotonic  decrease  in  the  apparent  rms  error  with 
depth  for  OTIS. 

Note  the  large  apparent  rms  errors  for  OTIS  between 
25  and  75  m  depth  in  the  midlatitude  eastern  Pacific 
compared  to  TEOTS  (  Area  1.  Fig.  10a).  From  ex¬ 
amination  of  the  apparent  mean  errors  for  this  area 
derived  from  the  same  bathy  dataset  (not  shown),  it 
is  apparent  that  the  OTIS  thermal  fields  were  biased 
cold  in  this  depth  range.  This  implies  that  the  OTIS 
MED  was  biased  shallow  in  the  region,  suggesting  room 
for  improvement  in  the  TOPS  mixed-layer  model  and/ 
or  the  manner  in  which  OTIS  analyzes  MED. 

Except  in  the  midlatitude  western  Atlantic  and  the 
Indian  Ocean,  where  the  apparent  errors  for  OTIS  in¬ 
dicate  the  largest  relative  improvement  over  TEOTS. 
the  verification  statistics  show  \ery  little  difference  be¬ 
tween  the  accuracy  of  the  two  products  at  the  surface. 
Improvements  of  OTIS  over  TEOTS  are  most  pro¬ 
nounced  below  about  100-150  m  depth.  This  is  con¬ 
sistent  with  the  conclusion  of  section  2  that  the  OTIS 
and  TEOTS  SST  fields  are  quite  similar,  w  hile  the  1 200 
and  1400  fields  exhibit  substantial  differences. 

e.  Depili-averat’ed  apparent  error  statislies 

Depth-averaged  apparent  rms  error  statistics,  cal¬ 
culated  by  integrating  the  rms  errors  of  Eig.  10  over 
various  depth  intervals  and  then  dividing  the  result  by 
the  depth  interval,  prov  ide  a  convenient  representation 
of  the  verification  results.  Depth-averaged  statistics  for 
I  ebruary -April  1988  are  presented  in  Fable  3  for  each 
of  the  nine  verification  areas.  Included  in  the  tables 
are  depth-averaged  apparent  rms  errors  for  OTIS  and 
TEOTS  and  the  percent  improvement  of  the  depth- 
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Fie  i .  10.  Apparent  root-mean-square  (  rms )  errors  for  ( )  I  IS  ( solid  line  I  anil  F  POTS  ( dashed  line  I  as  a  lu  net  ion 
of  depth  made  in  February -April  19X8.  The  error  bar  indicates  an  approximate  9tv;  confidence  interxal  lat  In 
the  midlalitude  eastern  Pacific  I  Area  I  )  from  97(1  independent  baths  reports,  t  b )  for  the  midlatitude  western  Pacific 
(  Area  2  )  and  1 2.32  hathys.  (  e )  for  the  midlatitude  eastern  Atlantic  (  Area  3  l  and  2  304  bathss.  ( d  I  for  the  midlalitude 
western  Atlantic  ( Area  4  )  and  2161  bathys.  I  e )  lor  the  Indian  Ocean  (  Area  3  l  and  .372  bathss.  ( I  )  for  the  tropical 
eastern  Pacific  ( Area  I A )  and  1 295  bathss.  ( g )  for  the  tropical  western  Pacific  I  Area  2  A 1  and  394  bathss.  t  h  l  for 
the  tropical  eastern  Atlantic  ( Area  3.4  )  and  104  bathss.  and  ( i )  for  the  tropical  western  Atlantic  (  Area  44  I  and 
948  bathys. 
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TEDTS 

OTIS 


TEOTS 

OTIS 


100 


X 

a. 

ID 

o 


200 


300 


RMS  ERROR  (DEC  C> 

1  2  3 


RMS  ERROR  (DEG  O 

1  2  3 

- 1 


CM  LIMIT 


RMS  ERROR  <0EG  C) 


Fl(i.  K).  [Conn maul) 


averaged  OTIS  rms  error  over  that  of  TEOTS  for  depth 
intervals  of  0-300  m,  0-150  m.  and  150-300  m. 

The  statistics  in  Table  3  show  regional  differences, 
with  the  apparent  rms  errors  for  both  OTIS  and  TEOTS 
largest  in  the  midlatilude  western  Atlantic  (Area  4) 
and  midlatitude  western  Pacific ( Area  2 ).  This  reflects 
the  influence  of  the  Gulf  Stream  and  Kuroshio  Current 
systems  and  their  related  mesoscale  eddy  fields.  The 
relative  improvement  of  OTIS  over  TEiOTS  is  largest 
in  the  Indian  Ocean  (  Area  5  ).  OTIS  performs  poorest 
relative  to  TEOTS  in  the  midlatilude  eastern  Pacific 
(Area  I  )  because  of  its  shallow  bias  in  MET)  noted 
above,  but  hath  OTIS  and  TEOTS  are  extremely  ac¬ 
curate  there  with  low  apparent  rms  errors  relative  to 
the  other  areas  and  climatology.  When  averaged  over 


all  areas  (9980  bathy  reports).  OTIS  shows  a  7'7  im¬ 
provement  over  TEOTS  in  0-300  m  apparent  rms  er¬ 
ror.  Stratified  by  depth  ranges.  OTIS  shows  only  a  3'7 
improvement  over  TEOTS  in  0-150  m  apparent  rms 
error,  but  an  I  I improvement  over  TEOTS  in  1 50— 
300  m  apparent  rms  error. 

It  is  useful  to  compare  the  accuracy  of  OTIS  to  that 
of  the  ocean  thermal  climatology  used  as  the  first-guess 
field.  Depth-averaged  apparent  rms  errors  for  OTIS 
and  this  climatology  for  February -April  in  the  nine 
verification  areas  are  shown  in  Table  4.  Also  included 
in  the  table  is  the  percent  improvement  of  OTIS  over 
the  climatology.  Averaged  over  all  areas.  OTIS  shows 
a  12'7  improvement  over  climatology  in  0-300  m  ap¬ 
parent  rms  error.  The  improvement  relative  to  cli- 
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Tabi  l  3.  Depth-averaged  apparent  rms  error  statistics.  [( n'  )  ')1’  \  for  OTIS  and  (  LOTS  for  February -April  1988. 


Area 

Area  definition 

Number  OT  IS  rms 

ofobs  error  (C) 

TEOTS  rms 
error  (C) 

Percent  improvement 
of  OTIS  over  TEOTS 

Depth  interval 

0-300  m 

1 

Midlatitude  LPAC 

97(1 

0.92 

0.88 

4 

*> 

Midlatitude  WPAC 

1232 

2.14 

2.17 

1 

.1 

Midlatiiude  ELAN  I 

2304 

1.14 

1.20 

5 

4 

Midlatitude  W1  ANT 

2101 

2.06 

2.40 

14 

s 

Indian  CXean 

.372 

1.17 

1.42 

18 

IA 

t  ropical  EPAC 

1295 

1.31 

1.37 

4 

2  A 

T  ropical  WPAC 

594 

1.35 

1.41 

4 

3  A 

Tropical  FLAN  1 

104 

1.12 

1.17 

4 

4  A 

Tropical  W1  AN  I 

948 

1 .08 

1.17 

8 

Average  over  all  areas 

1.36 

1.46 

7 

Depth  interval 

0-/50  m 

1 

Midlatitude  EPAC 

970 

1 .0.3 

0.92 

12 

A 

Midlatitude  WPAC 

1232 

2.1 1 

2.12 

1 

3 

Midlatitude  El  ANT 

2304 

1.08 

1.13 

4 

4 

Midlatitude  WI  AN  1 

2101 

2.24 

2.54 

12 

5 

Indian  Ocean 

372 

1.31 

1.48 

12 

1  A 

Tropical  LPAC 

1 295 

1.44 

1.44 

0 

2  A 

I  ropical  WPAC 

594 

1.34 

1.30 

3 

3  A 

Tropical  IT  AN  1 

104 

1.49 

1.44 

3 

4  A 

Tropical  W  L  AN  1 

948 

1.12 

1.18 

6 

Average  over  all  areas 

1.46 

1.51 

7 

Depth  interval 

130-300  m 

1 

Midlatitude  EPAC 

970 

0.80 

0.83 

4 

A 

Midlatitude  W  PAC 

1232 

2.17 

2.21 

*> 

3 

Midlatitude  ELAN  I 

2304 

1.19 

1.27 

7 

4 

Midlatitude  WI  AN  I 

2161 

1.87 

2.26 

18 

5 

Indian  Ocean 

372 

1.0.3 

1.36 

24 

LA 

I  ropical  EPAC 

1295 

1.19 

1.29 

8 

2  A 

Tropical  W  P  AC 

594 

1.37 

1.52 

It) 

3  A 

1  ropical  IT  AN  I 

104 

0.75 

0.90 

16 

4  A 

Tropical  WI  AN  1 

948 

1.04 

1.15 

10 

Average  over  all  areas 

1.27 

I  42 

1  1 

matology  is  largest  in  the  midlatitude  eastern  Pacific 
(Area  I  ).  and  smallest  in  the  tropical  western  Atlantic 
( Area  4 A  ). 

d  Depth-avcragi'd  grid-scale  error  statistics 

It  may  seem  surprising  that  the  differences  in  ap¬ 
parent  rms  error  between  O  f  IS.  I  LOTS,  and  clima¬ 
tology  are  not  larger,  given  the  qualitative  differences 
between  the  subsurface  thermal  fields  noted  in  section 
5.  To  explain  this,  we  write  expressions  for  the  apparent 
i  ms  vTi\*is  [  ( ./)" )" '  for  O' I  IS.  MO  IS,  and  climatol¬ 
ogy  as 

[  (  lolls)'  [(A'idiis  t  in’’)'  t  <it\mvl  (50) 

l  (  <r'  )  i  i  ois]  “  l  (  )i  km  s  1  (rr1)'  )  <it\m\l 

(31  ) 


[(<t‘)<T|.m]'  [(frU)(HM  +  (ff")'  +  <baihy)° 

(32) 

where  the  overbar  indicates  an  average  over  all  bathy 
data  in  a  verification  area,  and  the  grid-scale  error  a". 
the  subgrid-scale  error  a",  and  the  bathy  instrumental 
error  < h a  1 1 n  are  assumed  to  be  uncorrelated.  With  the 
thermal  fields  produced  by  the  two  models  and  cli¬ 
matology  all  defined  on  the  same  horizontal  grid,  dif¬ 
ferences  in  apparent  rms  error  [(a')  ]"'  can  be  due 
only  to  differences  in  a".  However,  for  the  grid  reso¬ 
lution  used  here  (about  320  km  ).  n1'  a"  s>  <„Ann  . 
Thus,  the  apparent  rms  error  statistics  are  highly  con¬ 
taminated  by  subgrid-scale  noise  ( <r")  which  tends  to 
overwhelm  and  mask  the  true  grid-scale  error  of  the 
fields  (n“)  and.  hence,  the  differences  in  the  relative 
performance  of  the  models. 
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Tabu  4 

Depth-averaged  apparent  rms  error  statistics.  [( «')•’ )"'.  for  O  l  IS  and 

elimatologv  for  1 

ebruary-  April  1988. 

Number  OTIS  rms 

CUM  rms 

Percent  improvement 

Area 

Area  definition 

of  obs 

error  (C) 

error  (C) 

of  OTIS  over  CUM 

Ih-plh  interval 

0-J00  m 

1 

Midlatitude  I  PAC 

970 

0.92 

1,24 

26 

2 

Midlatitude  WPAC 

1232 

2.14 

2.23 

4 

3 

Midlatitude  FI  ANT 

2304 

1.14 

1.21 

6 

4 

Midlatitude  WLANT 

2161 

2.06 

2.57 

20 

5 

Indian  Ocean 

372 

1.17 

1.22 

4 

1A 

Tropical  EPAC 

1295 

1.31 

1.57 

17 

2A 

Tropical  WPAC 

594 

1.35 

1.60 

16 

3A 

Tropica'  PLANT 

104 

1.12 

1.21 

7 

4A 

Tropical  WLANT 

948 

1.08 

1.09 

1 

Average  over  all  areas 

1.36 

1.55 

12 

Depth  interval 

0-150  m 

1 

Midlatitude  EPAC 

970 

1.03 

1.38 

26 

i 

Midlatitude  WPAC 

1232 

2.11 

2.14 

1 

3 

Midlatitude  FI.ANT 

2304 

1  08 

1.15 

5 

4 

Midlatitude  W1  ANT 

2161 

2.24 

2  65 

16 

5 

Indian  Ocean 

372 

1.31 

1.28 

*> 

1  A 

Tropical  EPAC 

1295 

1.44 

1.62 

1  I 

2  A 

Tropical  WPAC 

594 

1.34 

1.41 

s 

3  A 

Tropical  F!  A  N  I 

104 

1.49 

1.49 

0 

4  A 

Tropical  W!  AM 

948 

1.12 

1.04 

7 

Average  over  all  areas 

1.46 

IM 

7 

Depth  intend/ 

IM)-.iot>  in 

1 

Midlatitude  I  P  AC 

97() 

0.80 

1.10 

S7 

s 

Midlatitude  WP  AC 

1232 

2.17 

2.32 

6 

3 

Midlatitude  IT  AN  1 

2)04 

1.19 

1.27 

7 

4 

Midlamude  Wl  AN  1 

2161 

I.X7 

2  49 

25 

5 

Indian  Ocean 

372 

1.0) 

115 

1  1 

1  A 

I  Topical  1  PAC 

1 295 

1.19 

1.53 

s  s 

2  A 

I  ropical  WP  AC 

594 

1.37 

1.79 

24 

3  A 

1  ropical  II  AN  1 

104 

0."5 

0  94 

20 

4  A 

I  ropical  Wl  VN  1 

04  X 

104 

1.14 

9 

Average  over 

ill  areas 

1  27 

15) 

17 

Contamination  of  ocean  thermal  verification  statis¬ 
tics  with  observational  noise  is  a  fairly  general  problem, 
even  for  tiner-resolution  models.  \  technique  for  es¬ 
timating  the  uncontaminated  grid-scale  errors  of  a 
model  would  be  quite  useful,  and  can  be  derived  from 
the  formalism  presented  above,  f  rom  (  30 )  —  (  32  ). 

( 1  n  ‘‘  Ifnis]"  [  f  n"  )<  i  i\i  ■  (a' )o||S  (  n  him) 

(33) 

[  (  n“  I'll  ois]  [  (  n  ‘  )i  l  im 

t  ( n'  i  n  ms  ( n'  h  i  im  ]  (  34 ) 

The  quantity  (<i");  iim  is  the  mean-square  departure 
of  the  thermal  field  from  climatology  ( i.e..  the  mean- 
square  of  the  anomaly  held  1  calculated  over  the  loca¬ 
tions  of  the  bathy  observations  in  each  verification  area. 
It  can  be  estimated  from  the  analyzed  anomalies  during 


the  study  period.  Thus,  all  of  the  terms  on  the  right- 
hand-side  of  (  33  l  and  (34)  can  be  obtained  to  yield 
estimates  of  the  grid-scale  rms  errors  of  the  models 
( [  (  n" ) €*> i  is )'  '  and  [  (  n‘‘  In  ois]"  )  uneontammated  by 
observational  noise. 

Depth-averaged  values  of  the  grid-scale  rms  errors 
[  (  a“ ):  ]" '  lor  OHS  and  I IX)  IS  calculated  from  (  33  ) 
and  (  34 )  are  presented  in  I  able  5.  I  or  all  depth  inter¬ 
vals  and  areas,  the  grid-scale  rms  errors  of  fable  5  are 
smaller  than  the  apparent  rms  errors  of  fable  3.  re¬ 
flecting  the  removal  of  observational  noise.  Also,  the 
grid-scale  errors  do  not  as  clearly  take  on  relative  max¬ 
ima  in  the  western  boundary  current  areas  (  Areas  2 
and  4)  as  the  apparent  errors.  I  his  again  rellects  the 
suppression  of  observational  noise,  which  is  particularly 
strong  m  these  areas  due  to  the  presence  of  strong  and 
unresolved  mesoscale  eddies.  Comparison  of  the  grid- 
scale  errors  rather  than  the  apparent  errors  allows  a 
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Voi  t  mi  7 


1 

\Bl  l  5.  IXplliavcraged  grid-scale 

rms  error  statistics.  j(  a 

for  ( )  1 

ITS  and  HOTS  tor  I 

ehruarv  -April  10XX. 

Number  O  IIS  rms 

I  tOI  S  rms 

Percent  improvement 

Area 

Area  definition 

ofohs  c 

rror  (Cl 

error  ((  T 

off  ms  over  no  IS 

Ih'/’lh  iillcruil 

U-.HHi  m 

i 

Midlulilude  I  P  AC 

</7o 

0.55 

0  4X 

14 

t 

Midlatitude  WPAC 

1232 

1.00 

1.06 

6 

A 

Midlalitudc  l  I  \N  1 

2304 

(1.5S 

0.6l> 

16 

4 

Midlatitude  V\  1  AN  1 

2161 

(1.73 

1.43 

40 

s 

Indian  Ocean 

>7  1 

o.w 

1.06 

>  A 

1  A 

1  ropical  1  PA( ' 

12^5 

OX  1 

1)5)0 

10 

tropical  W  1’  \( 

3*14 

05)6 

1,06 

10 

y\ 

I  ropteal  1  1  AN  1 

1(14 

0.62 

03’  1 

13 

4  \ 

1  ropteal  \V1  ANT 

>)4S 

o  53 

0.60 

2  A 

A  v  erage 

over  all  areas 

(1.72 

0.00 

20 

Ih'/nli  inu  ruil 

0-/30  1)1 

1 

Midlatitude  I  I’  AC 

S)7(  I 

0  61 

0.30 

56 

a 

Midlatitude  VYPAC 

12  0 

1.12 

1.13 

5 

3 

Midlatitude  1  1  \\l 

2304 

0.20 

0.56 

44 

4 

Midlatitude  XVI  AN  1 

2161 

(1.03 

1.40 

5* 

5 

Indian  Ocean 

A  -*  "> 

0.S7 

l.l  1 

A  1 

1  A 

1  ropical  I  I’  AC 

1  2^5 

0.02 

0.02 

0 

:\ 

1  ropical  WPAC 

3*14 

1.02 

05>K 

4 

3  A 

1  ropteal  II  AN  1 

104 

O.S.S 

(1.75 

1  1 

4  A 

tropical  vs  1  AN  I 

‘MX 

0.62 

0.73 

13 

Average  over  all  areas 

0.70 

OSS 

10 

/Vpf/i  iun  nul 

{ '(l- in  m 

1 

Midlatitude  11’  AC 

*)7() 

0.4S 

0.33 

10 

A 

Midlatitude  WI’AC 

1232 

0.S5 

05)5 

1 1 

3 

Midlatitude  I  I  AN  I 

2304 

0.M0 

1.00 

1 1 

4 

Midlatitude  Wl  AN  1 

2IM 

0.44 

1.55 

6S 

5 

Indian  Ocean 

372 

0.4S 

t  .01 

33 

1  A 

1  ropical  1  P  AC 

1 203 

0.6S 

0.X3 

20 

2  A 

1  ropteal  VS 1’  At 

504 

0.03 

1.14 

IX 

3  A 

Tropical  IT  AN  1 

104 

0  43 

0,66 

32 

4  A 

1  ropical  VVT  ANI 

‘MX 

0.43 

0.66 

33 

Average  over  all  areas 

0.62 

0  00 

31 

more  meaningful  ranking  of  model  skill.  The  0-300  An  operational  test  of  OTIS  was  conducted  during 
m  grid-scale  rms  errors  of  Fable  5  averaged  over  all  F  ebruary.  March,  and  April  of  1988.  This  involved 
areas  show  O  F  IS  w  ith  a  20' V  improvement  over  running  OTIS  daily  in  parallel  w  ith  and  on  the  same 
TKOTS.  three  times  that  calculated  from  the  apparent  grid  as  the  TKOTS  analy  sis  ( the  existing  operational 
rms  error  statistics  of  Table  3.  product),  comparing  output  from  the  two  models 

qualitatively,  and  validating  the  output  of  the  models 
7.  Summary,  conclusions,  and  future  work  against  bathy  data  that  was  independent  ot  the  models. 

Qualitative  comparison  of  SS  I  fields  produced  by 
O  F  IS  is  an  ocean  thermal  analy  sis  system  designed  OTIS  and  TKOTS  shows  little  difference  between  the 
for  operational  use  at  f  NOC.  It  is  based  on  the  opti-  two  products.  Comparison  of  subsurface  fields,  how- 
mum  interpolation  data  assimilation  technique  and  ever,  shows  large  and  important  differences  between 
functions  in  an  analysis-prediction-analysis  data  as-  the  two.  OHS  subsurface  temperature  fields  exhibit 
similation  cycle  with  the  TOPS  mixed-layer  model,  sharper  and  more  realistic  horizontal  gradients,  more 
OTIS  provides  a  rigorous  framework  for  combining  physically  reasonable  large-scale  anomaly  patterns,  a 
real-time  data  ( e  g.,  ship,  buoy .  bathy.  and  satellite  ob-  better  ( though  still  very  poorly  resolved )  representation 
servalions).  climatology,  and  predictions  from  nu-  of  the  Gulf  Stream  and  Kuroshio  Current  systems,  and 
mcrical  ocean  prediction  models  ( eg..  FOPS )  to  pro-  a  better  representation  of  the  I  oop  Current  in  the  Ciulf 
duce  a  large-scale  sy  noptic  representation  of  ocean  of  Mexico. 

thermal  structure.  Validation  of  OTIS  and  TfX)  IS  against  independent 
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bathy  data  (i.c..  unassimilated  into  the  analyses  at  the 
time  of  validation)  indicates  that  OTIS  produces  a 
more  accurate  representation  of  ocean  thermal  struc¬ 
ture  than  TEOTS.  Apparent  rms  errors  for  OTIS  arc 
generally  less  than  those  for  TEOTS,  particularly  below 
about  100-150  m  depth.  The  0-300  m  depth-averaged 
apparent  rms  errors,  calculated  from  a  90-day  accu¬ 
mulation  of  data  and  averaged  overall  areas  to  produce 
a  sample  size  of  9980  bathy  reports,  show  OTIS  with 
a  7%  improvement  over  TEOTS.  The  0-300  m  depth- 
averaged  grid-scale  rms  errors,  estimated  via  a  tech¬ 
nique  to  remove  the  effect  of  observational  noise  from 
the  statistics  and  averaged  over  all  areas,  show  OTIS 
with  a  20%  improvement  over  TEOTS. 

OTIS  was  designated  as  the  Navy’s  new  operational 
global-scale  ocean  thermal  analysis  in  July  1988.  Pres¬ 
ent  and  future  work  involves  the  application  of  a  new 
and  more  sophisticated  version  of  OTIS  on  eddy-re¬ 
solving  regional  grids,  and  the  continued  monitoring 
of  verification  statistics  derived  from  independent  bathy 
data. 
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APPENDIX  A 

Definition  of  Symbols 

Ak  Inverse  of  east-west  decorrelation  scale  of 

resolvable  thermal  anomalies  at  gridpoint 
k. 

Bk  Inverse  of  north-south  decorrelation  scale 

of  resolvable  thermal  anomalies  at  grid- 
point  k. 

C\  Inverse  of  temporal  decorrelation  scale  of 

resolvable  thermal  anomalies  at  gridpoint 
k. 

e  Mean  of  instrumental  error  (i.c.,  bias)  for 

MCSST  observations. 

(I a  Rms  TOPS  error  growth  over  one  forecast 

cycle  at  gridpoint  k. 

Iik"  Analyzed  mixed-layer  depth  at  location  /. 

Ilf  Climatological  mixed-layer  depth  at  grid- 
point  k. 

h,"  Observed  mixed-layer  depth  at  gridpoint  k. 

M  Number  of  MCSST  observations  uveraged 

around  a  gridpoint  to  form  a  super-ob¬ 
servation. 

N  Number  of  observations  assimilated  at  grid- 

point  k. 

T  Temperature. 

T’  Reference  climatological  SST  gradient. 

Tk"  Analyzed  temperature  at  gridpoint  k. 


Tic c  Climatological  temperature  at  gridpoint  k. 

T,"  Observed  temperature  at  location  /. 

Tkp  TOPS  predicted  temperature  at  gridpoint  k. 

-  Depth,  measured  positive  downward  from 

the  surface  in  meters. 

Zn  Deepest  reported  depth  in  a  bathy  obser¬ 

vation. 

Tj  Individual  MCSST  observation  at  location 

/. 

<u,  Weight  given  to  observation  /  assimilated  at 
gridpoint  k. 

fik  Weight  given  to  TOPS  prediction  at  grid- 

point  k. 

yk  Weight  given  to  climatology  at  gridpoint  k. 

j)i ,  Correlation  between  resolvable  thermal 

anomalies  at  location  /  and  location  j  in 
space-time, 

X,"  Noise-to-signal  ratio  for  observation  /. 

Xf  Noisc-to-signal  ratio  for  TOPS  prediction  at 
gridpoint  k. 

X"  Noise-to-signal  ratio  for  subgrid-scale  error 

at  location 

X;:„  Average  of  X/’  over  entire  horizontal  analysis 
domain. 

A  A'  Reference  grid  space. 

A.v„  East-west  distance  between  location  /  and 

location  j. 

Ay  a  North-south  distance  between  location  i  and 
location  j. 

Ala  Time  difference  between  observation  /  and 
observation  j. 

(A7V’)«  Change  in  TOPS  predicted  temperature 
from  level  zz  —  1  to  level  n  (level  0  implies 
the  surface)  at  gridpoint  k. 

of  Rms  observational  error  of  observation  /. 

Okp  Rms  error  of  TOPS  prediction  at  gridpoint 

/  V  . 

of  Long-term  rms  variation  of  resolvable  ther¬ 
mal  anomalies  about  climatology  at  lo¬ 
cation  i. 

of  Rms  difference  between  bathy  observations 
and  climatology  at  location  /. 

rr*t,  Rms  difference  between  bathy  observations 
and  climatology  calculated  over  entire 
horizontal  analysis  domain. 

of  Expected  rms  temperature  difference  be¬ 
tween  Tf  and  independent  bathy  obser¬ 
vations. 

of  Rms  subgrid-scale  error  at  location  i. 

of  Rms  error  added  to  MCSST  super-obser¬ 

vation  as  a  result  of  the  way  in  which  the 
block  averaging  of  the  MCSST  data  is 
clone. 

of  Rms  error  of  resolvable  thermal  field  in  re¬ 
sulting  analysis  at  gridpoint  k. 

(  ^initial  Rms  error  of  the  resolvable  thermal  field 
from  the  previous  analysis  initializing 
TOPS  at  gridpoint  k. 
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(<t/’)ex  Rms  error  added  to  / th  bathy  observation 
below  depth  :B  by  the  downward  extrap¬ 
olation  of  the  bathy. 

£  Standard  deviation  of  instrumental  error  for 

MCSST.  SHIP,  or  BATHY  observations. 
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